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Abstract. We review the observational evidences and the possible theoretical expla- 
nations for the abundance gradients in the Galactic disk. In particular, we discuss the 
implications of abundance gradients and gradients of abundance ratios on the mecha- 
nism for the formation of the Galaxy. We conclude that an inside-out formation of the 
Galaxy, and in particular of the Galactic disk, where the innermost regions are assumed 
to have formed much faster than the outermost ones, represents the most likely expla- 
nation for abundance gradients and we predict that the abundance gradients along the 
Galactic disk have increased with time. 



1 Introduction 

The existence of abundance gradients along the Galactic disk seems now to be 
well established with agreement between different data sources. 

— HII regions: Shaver et al. (1983) first measured abundance gradients from 
Galactic HII regions indicating that O/H, N/H, Ar/H show gradients of the 
order of —0.07 dex Kpc -1 between 4 and 14 Kpc along the Galactic disk. 
More recently, Vilchez and Esteban (1996) found a flattening of O/H, N/H 
and S/H gradients for Rqc > 12 Kpc suggesting that these gradients present 
a bimodal behaviour. 

— Planetary Nebulae: Maciel and Chiappini (1994) and Maciel and Quireza 
(1998) found O/H, Si/H, Ne/H, Ar/H gradients from planetary nebulae very 

similar to those found from HII regions, namely between —0.04 0.07 dex 

Kpc -1 in the galactocentric distance range 4-14 Kpc. 

— B type stars: Kaufer et al. (1994) found that B type stars did not exhibit 
any appreciable gradient but more recent studies (Smartt and Rolleston, 
1997 and Gummersbach et al. 1998) found a gradient of oxygen of —0.07 ± 
0.01 dex Kpc -1 between 6 and 18 Kpc and — 0.07± 0.02 dex Kpc -1 between 
5 and 14 Kpc, respectively. These results agree with the gradient estimated 
by Shaver et al. (1983) from HII regions and with the planetary nebula data. 

— Open clusters: as indicated in the review by Friel (1995) the iron abundance 
in open clusters shows a gradient of [Fe/H] of — 0.095± 0.017 dex Kpc -1 
between 7 and 16 Kpc but Friel (1998) showed a more recent determination 
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of the [Fe/H] gradient shallower than the previous one and closer to the O 
gradients as estimated recently from B type stars . 

It is therefore encouraging that different data sources seem to agree at least 
on the value of the oxygen gradient. However, we would be careful in concluding 
that the gradients of all the other heavy elements (C, N, Fe, Ar. etc..) should 
be the same as oxygen since from a theoretical point of view we would expect 
small differences due to the different nucleosynthetic production mechanisms and 
different progenitors of the elements, as we will see in the next paragraphs. 

The question we should ask now is : how can we obtain abundance gradients 
along the Galactic disk? 

Several mechanisms have been proposed in the last years: 

— i) Variations in the gas fraction as a function of the galactocentric distance 
Rgc (original idea of the "simple model") but the predicted gradients are 
too flat compared to observations (Pagel 1989). 

— ii) Variations in the initial mass function (IMF) but most of the suggested 
variations (massive stars increase with decreasing metallicity Z) tend to de- 
stroy the gradient! Carigi (1996) tried the opposite case, where low mass 
stars increase with decreasing metallicity, and showed that steep abundance 
gradients can be obtained but that the same model does not fit the chemical 
evolution of the solar ncighbourhhod. However, before concluding that there 
is no viable solution with a variable IMF we should try to test the variation 
of the IMF recently proposed by Larson (1998) who suggested a variation 
in the lower mass cut-off with galactocentric distance rather than a variable 
slope. 

— iii) Metal dependent stellar yields a la Maeder (1992) (less metals are pro- 
duced at higher Z) produce shallow gradients along the Galactic disk, as it 
has been shown by several papers (Giovagnoli and Tosi, 1995; Carigi 1996), 
which do not agree with the most recent observational estimates. 

— iv) Variations in the efficiency of star formation (SFR = v(Rcc)o~g a s)- 
Prantzos and Aubert (1995) explored this possibility by assuming v cx Rqq 
and found shallow gradients, again too shallow if compared with the most 
recent data. 

— v) Variations of the star formation rate (SFR) relative to the infall rate (IR) , 
namely SFR/IR, with Rgc- This hypothesis is also known as biased infall 
since one can obtain such a situation by assuming that the Galactic disk 
formed by infall of gas occurring at a faster rate in the innermost regions 
than in the outermost ones. The physical reason for a biased infall can be 
found in the fact that the gas tends to collapse more quickly in the center of 
the spheroid so that a gradient in the gas density, being lower in the more 
external parts, is established. In this situation the gas which continues to 
infall towards the disk falls more efficiently towards the center than towards 
the external regions, due to the stronger tidal force acting at the center (see 
also Larson 1976). Several authors showed that the biased infall can well 
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reproduce steep abundance gradients along the disk, especially if coupled 
with a star formation rate proportional to some power k of the gas density 
SFR oc <r k gas with k > 1 (Tosi 1988; Matteucci and Francois 1989; Matteucci 
et al. 1989; Chiappini et al. 1997). Therefore, we consider this as a viable 
solution to explain the existence of abundance gradients. 

— vi) Radial flows along the Galactic disk arc not enough by themselves to 
produce steep gradients but they work well if coupled with biased infall or 
a threshold in the SFR under specific conditions (Clarke 1989; Sommer- 
Larsen and Yoshii 1989; Goetz and Koeppen 1992; Chamcham and Tayler 
1994, Edmunds and Greenhow 1995). 

— vii) A biased outflow assuming a stronger and earlier wind in the external 
regions than in the internal ones can also produce steep abundance gradients, 
as shown for elliptical galaxies by Martinelli et al. (1998). At this point one 
could ask : how can we distinguish between biased infall and biased outflow?. 
There is a clear way to do it: thru the abundance ratios. In fact, under 
the biased infall hypothesis the predicted present time [a/Fe] ratios tend to 
decrease with galactocentric distance. This is due to the fact, as shown by 
Matteucci (1992), that if a region of the Galaxy evolves slower than another 
one then, at the same [Fc/H], the slow region will show a lower [a/Fe] ratio, 
due to the longer timescale required to achieve the same [Fe/H] relative to the 
fast region and therefore to the larger contribution from type la supernovae 
which are responsible for the bulk of iron production (Greggio and Renzini, 
1983; Matteucci and Greggio 1986), at that given [Fe/H]. On the other hand, 
under the biased outflow hypothesis Martinelli et al. (1998) predicted that 
the [a/Fe] ratio increases with the galactocentric distance, due to the fact 
that the galactic wind occurs before in the external regions than the internal 
ones, and that it is likely that no more or negligible star formation has 
taken place after the wind event. As a consequence of this, the stars in the 
external regions will not show the iron contributed by supernovae of type la 
but rather that produced by massive stars. In order to decide which one is 
the mechanism responsible for the abundance gradients in the Galactic disk 
one should check the gradient of the [a/Fe] ratio. 

In the following sections we will describe examples of biased infall models 
for the chemical evolution of the Galactic disk and compare their predictions 
with the available data. In particular, we will discuss abundance gradients of 
several chemical elements and of their ratios and by comparing these results 
with observations we will try to infer some conclusions on the mechanism of 
formation of the Galaxy. 

2 Chemical Evolution Models 

The Matteucci and Francois (1989) (hereafter MF) model assumes that the Galac- 
tic disk is divided in several independent rings 2 Kpc wide with no exchange of 
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matter between them. It assumes also that the disk forms out of infall of pri- 
mordial gas which has accumulated with different timescales at different galac- 
tocentric distances. The galactic halo represents the earliest phase of this infall 
process and no clear distinction between halo and disk is present. In particular, 
the assumed law for the disk formation is a linear relation between r and Rgc'- 



It is worth mentioning that a quadratic or exponential law for r does not 
produce relevant differences, as shown by Carigi (1996), in the predicted abun- 
dance gradients. The adopted stellar yields arc independent of metallicity and 
the star formation rate is proportional to both the surface gas density and the 
total surface mass density, SFR oc o^ as u^ ot . Matteucci and Francois concluded 
that a good fit to the observed abundance gradients (Shaver et al. 1983) and gas 
distribution in the disk requires k=l.l and h=0.1. 

Chiappini's et al. (1997) (hereafter CMG) model is a two-infall model, in the 
sense that assumes that the halo and the thick disk formed much more quickly 
than the thin disk during a first infall episode, whereas the thin disk formed by 
means of a second completely independent episode of primordial infalling gas. 
In particular, the timescale for the formation of the halo-thick disk is th ~ 1 
Gyr and does not vary significantly with the galactocentric distance whereas the 
timescale for the thin disk (td) varies as: 



In this approach the evolution of the halo-thick disk is completely disentan- 
gled from that of the thin disk. The star formation rate is the same as in MF 
except that k = 1.5 and h = 0.5. Is is worth noting that the higher k value 
adopted by CMG together with the longer timescale for the thin disk formation 
(8 Gyrs for the solar vicinity) ensures a good fit of the reviewed G-dwarf metal- 
licity distribution by Rocha-Pinto and Maciel (1996). Moreover, this choice for 
the star formation exponent is in agreement with recent estimates of the star for- 
mation rate in spirals by Kennicutt (1998). This model also adopts a threshold 
surface gas density (a gas = 7M Q pc~ 2 ) for the star formation rate (Kennicutt, 
1989) which naturally produces a iatus in the star formation rate between the 
end of the thick disk phase and the beginning of the thin disk phase. Such 
a iatus seems to be real since it is observed either in the plot of [Fc/O] versus 
[O/H] (CMG) or in the plot of [Fe/Mg] versus [Mg/H] (Bernkopf and Fuhrmann, 
1998). As we will see, the assumed threshold is also responsible for the flatter 
than MF's gradients in the inner thin-disk regions and the very flat gradients 
in the outermost thin-disk regions. Chamcham and Tayler (1994) also found a 
bimodal gradient due to the adoption of a threshold in the star formation rate. 

2.1 Basic Equations 



t(R gc ) = QAQARgc - 1-59 



(1) 



td(Rgc) = 0.875R gc - 0.75 



(2) 



The basic equation in both models (MF and CMG) are: 
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In CMG the infall term is denned as: 

d -^U = (Xi)inf e -*/« + (X i)inf e-(*-«-)/™ (4) 

dt a(r,tG) cr(r,tG) 

where Gj(r, t)i n f is the normalized surface gas density of the infalling material in 
the form of the element i, {Xi) ini gives the composition of the infalling gas which 
is assumed to be primordial. t max is the time of maximum gas accretion onto the 
disk, and tjj and tjj are the timescales for the mass accretion in the halo and 
disk components, respectively. These are the two really free parameters of the 
model and are constrained mainly by comparison with the observed metallicity 
distributions for halo and disk stars. The t ma x value is chosen to be 2 Gyrs and 
roughly corresponds to the end of the halo-thick disk phase. The quantities A(r) 
and B(r) are derived by the condition of reproducing the current total surface 
mass density distribution in the halo and the disk, respectively. The current total 
mass distribution is taken from Rana (1991). 

In the previous model of MF the infall term was simply: 

dG t (r,t) ini A(r) t/r 
dt a{r, t G ) 

with a unique r for the halo and the disk and where A(r) was chosen under 
the condition of reproducing the present time total surface mass density in the 
disk. In both the models the adopted IMF is taken from Scalo (1986). 



3 Model Predictions 

A good model for the chemical evolution for our Galaxy should honour a minimal 
number of observational constraints both in the solar neighbourhood and in 
the whole disk. Both the models described above were tested to reproduce the 
following constraints in the solar vicinity: i) the relative number of disk and halo 
stars in the solar cylinder, ii) the present day gas fraction, iii) the type la and II 
supernova rates and their ratio, iv) the solar abundances, v) the age-metallicity 
relation, vi) the present day infall rate, vi) the G-dwarf metallicity distribution, 
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vii) the [el/Fe] versus [Fe/H] relations. Under these conditions MF concluded 
that good abundance gradients in agreement with the Shaver et al. (1983) data, 
which agree with the most recent estimates, are reproduced by assuming an 
inside out formation of the disk, as described by equation (1). In addition, 
they showed that the star formation law should have < k < 1.5 with a best 
value of k — 1.1 in order to reproduce, at the same time, abundance gradients 
and an acceptable distribution of the surface gas density along the disk. In 
fact, the gas distribution along the disk strongly depends on the assumed power 
of the surface gas density in the star formation law. The predicted gradients, 
in the galactocentric distance range 4-14 Kpc, are: ^a^J*^ = —0.065 dex 

Kpc- 1 ; Alo9 £^' H) = -0.066 dex Kpc" 1 ; Alog £f c /H) = -0.085 dex Kpc" 1 ; 

Al °9(Fe/H) = _ Q 0? ^ K 1 AlogCHe/H) = _ Q QQgg dcx R -1 Ag we can 
ZIKgc AHqc 

see, most of these gradients are in agreement with the most recent estimates. An 
interesting fact is that gradients of different elements are slightly different and 
this is due to the different nucleosynthesis history of each element. For example 
14 N, which is mostly produced in low and intermediate mass stars as a partly 
primary and partly secondary (proportional to the original abundances of O and 
C) element, shows a steeper gradient than oxygen which is mostly produced in 
massive short lived stars. In addition, the fact of having a secondary nature plays 
a role in defining the absolute value of the abundance gradient. The gradient of 
Fe is also slightly steeper than that of oxygen since iron is mostly produced on 
long timescales by type la SNe (white dwarfs in binary systems). Therefore, both 
the timescales and the primary /secondary nature of a chemical element play a 
role in defining its gradient along the Galactic disk. 

CMG predicted the existence of bimodal gradients along the Galactic disk, in 
the sense that their gradients are flatter in the outermost regions and steeper in 
the internal ones (R < 10 Kpc), in agreement with Vilchez and Esteban (1996). 
In particular, their gradients, in the range 4-14 Kpc, are: 



^ArJc^ = -°- 032 dcx K P C_1 in tllC innCr rC S ion and ^ArJc^ = -°- 018 dex 

Kpc -1 in the outer region, M ° 3 ^ lR ^ ( { H " > = —0.037 dex Kpc -1 in the inner region 

and Alo9 /^J H) = -0.028 dex Kpc" 1 in the outer region, Al °l^ H) = -0.04 

dex Kpc -1 in the inner region and Al °^^ H ' > = —0.017 dex Kpc -1 in the outer 
region. 

The gradients are all flatter than those predicted by MF and the main reason 
resides in the fact that CMG adopted a threshold surface gas density below 
which the star formation is assumed to stop. This mechanism acts mostly in the 
external regions where the amount of gas is almost always close to the threshold. 
However, the best model of CMG predicts a gas distribution along the disk (see 
figure 2) in better agreement than that predicted by MF and this is also due to 
the assumed threshold in the star formation rate. 
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An alternative explanation for the formation of abundance gradients is a variable 
IMF, as discussed in the introduction. In hgure 1 we show the predicted gradient 
of oxygen by the best model of CMG (Model A), by a model considering an IMF 
which predicts more massive stars in lower density regions (Padoan et al. 1997) 
(Model B) and the predictions of Model C where both the variable IMF and the 
variable star formation efficiency are considered. The star formation efficiency 
is assumed to increase with decreasing Rgc as in Prantzos and Aubert (1995). 
It is clear that a variable IMF of this kind tends to destroy the gradient unless 
the additional hypothesis of the variable star formation efficiency is considered. 
However, as shown in figure 2, this last model (Model C) does not fit at all the 
observed gas distribution along the disk and therefore it should be rejected. 
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Fig. 1. The predicted oxygen abundance along the Galactic disk by three different 
models 
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Fig. 2. Observed and predicted gas distribution along the Galactic disk. The data and 
the models are indicated in the figure. 



4 Discussion and conclusions 

In this paper we have shown that the most realistic scenario to explain abundance 
gradients along the Galactic disk is the inside-out formation of the disk, in the 
sense that the internal parts of the disk are assumed to have formed on shorter 
timescales than the external parts. This allows us to reproduce gradients in 
agreement with the observations as long as the exponent of the star formation 
law, k, is larger than unity. At the same time we can reproduce the observed gas 
distribution along the disk as long as k < 1.5. Therefore, we can envisage a range 
of possible values for k, namely 1.1 < k < 1.5. Better agreement is obtained if 
also a dependence on the total surface mass density is considered in the star 
formation law (see also Tosi 1988). 

Concerning the gas distribution along the Galactic disk, we showed that the 
assumption of a threshold in the star formation rate improves the agreement 
with observations although it predicts slightly too shallow gradients, at least in 
the light of the more recent data. 

The assumption of a variable IMF where massive stars are more abundant in 
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low metallicity, low density regions tend to destroy the gradient along the disk, 
at variance with observations, although we can not exclude the existence of a 
particular variation of the IMF which can reproduce all the main observational 
constraints. 

Therefore, our main conclusion is that an inside- out scenario with constant 
IMF is the most favoured for the formation of abundance gradients. 

In order to check this point we should try to observe [a/Fe] ratios in stars 
in the outermost disk regions. In fact, we have also shown that one should 
expect differences (although small) between the gradients of elements produced 
on different timescales. For this reason we expect that the [a/Fe] ratio has a 
negative gradient with increasing Roc, and future measurements of such a 
ratio at large galactocentric distances can assess this point. Recently, Nissen 
and Schuster (1997) identified metal rich halo stars ([Fe/H] ~ -1 dcx) with both 
high and low [a/Fe] ratios. The halo stars with low [a/Fe] ratios tend to be on 
orbits biased to the outer halo. This could be the consequence of the inside-out 
Galaxy formation, in the sense that the outer halo formed more slowly than the 
inner one in analogy with what happens in the disk. In particular, we predict 
that the stars in the outer disk should show this effect even more than the outer 
halo stars. 

Another important fact, strictly related to the mechanism of Galaxy forma- 
tion, is the evolution of gradients with time. On this point there is no general 
agreement between different data sources and different authors! 

From the observational point of view the situation can be summarized as 
follows: 

— Planetary nebulae seem to suggest that gradients steepen with time (Ma- 
ciel and Koeppen 1994); in fact, type III PNe (the oldest) have a flatter 
gradient than type II and I PNe, but dynamical effects could perhaps be 
responsible for this behaviour. 

— Open Clusters suggest no clear variation of gradients with time, and may 
be a flattening (Carraro et al. 1998), but again dynamical effects could have 
played a role. 

From the theoretical point of view some models predict a steepening of the 
gradients with time (Tosi 1988; MF; Goetz and Koeppen 1992; CMG) whereas 
others predict the opposite trend (Ferrini et al. 1994; Allen et al. 1998). Gener- 
ally, models with biased infall (inside-out formation of the disk), such as those 
discussed here, predict a steepening of gradients in time due mainly to the as- 
sumed long infall timescales at large galactocentric distances (see also the review 
of Tosi 1996). Clearly the evolution of gradients with time is strictly related to 
the mechanism of Galaxy formation and is an important point to assess. How- 
ever, given the situation, it is still premature to draw firm conclusions and more 
precise observational data are necessary. 
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